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ABSTRACT 


A  mathematical  model  of  a  small  solid  propellant  rocket  has 
been  programmed  for  use  with  a  Univac  Solid  State  90  Digital  Computer 
by  means  of  Fortran  I  and  Fortran  II.  In  addition  to  calculating 
the  pressure-time  and  thrust-time  transients,  the  program  computes 
time-averaged  pressure,  thrust  coefficient,  thrust,  flow  rate,  and 
burning  rate,  as  well  as  total  impulse  and  specific  impulse. 

The  report  includes  the  derivation  of  equations,  the  method  of 
solution,  the  program  in  Fortran  I  and  Fortran  II,  an  example  solu¬ 
tion  compared  with  experiment,  and  a  recommended  procedure  for  use 
of  the  computer  program  in  the  design  and  development  of  a  small 
solid  propellant  rocket  motor. 
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INTRODUCTION 


Basic  rocket  equations  and  empirically  derived  propellant  con¬ 
stants  and  motor  dimensions  may  be  used  to  predict  the  performance  of 
small  solid  propellant  rocket  motors.  Programming  these  equations  on 
a  digital  computer  reduces  the  solution  time  for  such  performance  es¬ 
timates  to  a  matter  of  minutes.  This  report  is  concerned  with  the 
preparation  of  a  suitable  digital  computer  pro'gram  for  this  purpose. 

Specifically,  this  analysis  is  intended  to  predict  the  pressure 
and  thrust  vs  time  transients  in  rocket  motors  that  are  large  enough 
to  assume  equilibrium  between  the  gas  generation  rate  and  the  gas 
discharge  rate,  and  small  enough  so  that  temperature  gradients 
in  the  propellant  grains  may  be  considered  negligible. 

The  objective  of  this  report  is  to  present  a  digital  computer 
program  which  will  serve  as  a  satisfactory  tool  in  the  design  and 
development  of  certain  small  solid  propellant  rocket  motors. 


ASSUMPTIONS 


The  following  assumptions  have  been  made  in  the  construction 
of  the  mathematical  model  of  the  solid  propellant  rocket  motor. 

1.  Initially  (time  =  0),  the  igniter  system  has  pressurized 
the  chamber  to  a  pressure,  P0. 

2.  All  propellant  surfaces  are  ignited  isochronically ,  and 
burning  is  directed  normal  to  all  exposed  propellant  surfaces. 

3.  Propellant  geometry  can  be  described  by  a  form  function, 
and  propellant  is  completely  burned  upon  disappearance  of  the  web. 

4.  The  space-averaged  temperature  of  the  combustion  products 
remains  constant  at  the  isobaric  (i.e.,  constant  volume,  adiabatic) 
flame  temperature,  T. 

5.  The  gas  generation  rate  and  the  gas  discharge  rate  are 
equal. 

6.  Temperature  gradients  in  the  propellant  grain  are  negli¬ 
gible. 


★See  REFERENCES. 
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EQUATIONS 


The  equations  predicting  the  pressure*  vs  time  relationship  are 
combined  in  Appendix  A  into  the  following: 


The  £hrust  value  corresponding  to  a  given  pressure  value  is  de¬ 
fined  by: 


F  =  <p  6  Cf  At  P  (Ref  2) 


in  which' 


and 


_2±L_ 


(Ref  3) 


Equations  1  to  4  are  used  to  predict  the  pressure  and  thrust  vs 
time  curves.  In  addition  to  predicting  these  relationships,  it  is 
desirable  to  compute  additional  parameters  that  characterize  a  spe¬ 
cific  solid  propellant  rocket  motor.  Additional  parameters  computed 
in  this  program  are : 


A.  Time-Averaged  Pressure 


♦Unless  otherwise  stated,  the  word  "pressure"  will  refer  to  head 
pressure . 
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B.  Time -Averaged  Thrust  Coefficient 


Cf  -  Of  -  Y  (6) 

C.  Time -Averaged  Thrust 

F  -  <P  6  Cf  P  At  (7) 

D.  Time -Averaged  Gas  Flow  Rate  Through  the  Nozzle 


cex  a  cd  At  p 

E.  Specific  Impulse 

lap  - 

F.  Total  Impulse 

^•tot  =  ^  *-c 

G.  Time -Averaged  Burning  Rate 

r  =  B  P n 


(8) 


(9) 


(10) 


(U) 


These  additional  parameters  are  computed  after  the  computation 
of  the  pressure  and  thrust  vs  time  relationships. 


METHOD  OF  SOLUTION 


Equation  1  is  an  implicit  expression  of  pressure,  P,  as  a  func¬ 
tion  of  the  fraction  of  web  burned,  z.  All  other  values  in  this  e- 
quation  are  considered  constants.  Since  z  ranges  from  0  to  1,  it  is 
conveniently  used  as  the  independent  variable.  For  a  given  value  of 
z  and  an  assumed  value  of  P,  the  right  hand  side  of  Equation  1  may  be 
calculated.  If  the  result  of  this  calculation  is  in  close  agreement 
with  the  assumed  value  of  P,  then  the  assumed  value  is  adequate.  If 
there  is  poor  agreement,  then  another  estimate  for  P  is  made  and  the 
calculation  is  repeated.  This  process  continues  until  good  agreement 
between  the  assumed  value  and  the  calculated  value  is  obtained.  If 
systematic  estimates  are  made,  the  number  of  iterations  are  reduced. 
The  specific  details  of  this  iteration  as  used  in  the  computer  pro¬ 
gram  are  found  in  Appendix  B . 
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For  each  value  of  P  computed  in  this  fashion,  a  corresponding 
value  of  F  is  calculated  by  means  of  Equation  2.  Before  Equation  2 
is  solved,  the  parametric  Equations  3  and  4  must  be  solved  for  Cf. 
Details  of  the  Cf  calculation  by  use  of  the  Newton-Rapthson  Method 
are  found  in  Appendix  C. 

An  explanation  of  the  <p  6  quantity  is  in  Appendix  D. 


COMPUTER  PROGRAM 


The  method  of  solution  described  above  and  in  Appendices  B  and 
C  is  the  basis  for  the  computer  program  which  was  written  in  Fortran 
I  and  Fortran  II  for  use  with  a  Univac  Solid  State  90  Computer.  These 
programs  appear  in  Appendix  E  of  this  report. 

*  To  test  this  computer  model,  parameters  describing  the  XM15 
Catapult,  an  experimental  solid  propellant  rocket  motor  designed  to 
eject  the  nose  capsule  of  a  fighter  aircraft,  were  used  as  input  data 
with  the  program.  These  parameters,  together  with  the  computer  re¬ 
sults  compared  to  experimental  results,  are  shown  in  Appendix  F. 


CONCLUSIONS 


1.  The  method  presented  is  useful  in  predicting  the  performance 
of  small  solid  propellant  rocket  motors  as  used  in  Propellant  Actuated 
Devices  (PAD)  . 

2.  Fortran  has  proved  effective  in  accomplishing  a  digital  com¬ 
puter  simulation  of  small  solid  propellant  rocket  motors. 


RECOMMENDED  PROCEDURE 


It  is  recommended  that  the  computer  program  presented  here  be 
incorporated  in  the  performance  estimation  and  early  design  and  de¬ 
velopment  phases  of  small  propellant  rocket  motors  as  used  in  PAD. 
The  procedure  is  as  follows: 

l.  Select  initial  design  parameters  consistent  with  envelope 
requirements . 
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2.  Perform  exploratory  computer  run  to  evaluate  initial  design 
parameters.  — -  - 

3.  Adjust  parameters  and  make  computer  tests  until  agreement 
between  the  computer  output  and  the  design  requirements  is  obtained. 

4.  Build  and  test  and  experimental  model  based  on  the  parameters 
selected  in  3  above. 

5.  Adjust  propellant  parameters  to  match  the  experimental  per¬ 
formance.  This  is  the  calibrated  analytical  model. 

6.  Use  calibrated  analytical  model  to  establish  optimum  design 
parameters . 


FUTURE  WORK 


Future  work  will  be  directed  toward  the 
for  the  automated  design  of  more  complex  PAD 


preparation  of  programs 
systems . 
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APPENDIX  A 


DERIVATION  OF  EQUATION  1 


Each  time  the  propellant  web  regresses  by  an  increment  dz  =  ilH, 
the  weight  of  charge  burnt  (cb)  ,  in  pounds,  is  given  by  w 


cb  «  z(l  +  0f)  cp 


(A-l) 


where  z  is  the  total  fraction  of  web  burnt  (dimensionless); 

f  is  the  total  fraction  of  web  remaining,  i.e.>  1  -  z 
(dimensionless) ; 

0  is  the  form  factor  (dimensionless) 

Cp  is  the  total  charge  weight  at  the  start  (lb)  . 

The  free  volume  (V),  in  cu  in.,  is  the  sum  of  the  initial  free 
volume  (VQ) ,  in  cu  in.,  and  the  volume  left  by  the  burned  propellant. 


V  =  vo  4*  — 

P 


(A-2) 


where  p  is  the  gravimetric  density  of  the  propellant  (lb/in.3) . 


The  rate  of  web  regression  (r)  in  inches  per  second  is  given  by 


r  -  BP*V 


(A-3) 


where  B  is  the  regression  rate  coefficient  (in./sec) (in.^/lb)n; 

n  is  the  regression  rate  exponent  (dimensionless) ; 

Pav  is  the  pressure  (lb/in. 2)  averaged  between  the  pressures  cor¬ 
responding  t»  k  and  z  +  dz,  i.e.,  Pav  ■  p(ztdz) 

The  time  (dt) ,  in  seconds,  required  for  the  web  to  burn  a 
fraction  dz  is  given  by 


dt  ■ 


w  dz 


(A-4) 
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where  w  is  the  total  web  (in.) . 


The  weight  of  gas  (cex) ,  in  pounds,  that  has  been  expelled  up 
to  time  t  is  given  by 


cex  ®  cd 


P  dt 


where  C^  is  the  discharge  coefficient  (sec“^); 

At  is  the  throat  area  of  the  nozzle  (in.2); 
P  is  the  pressure  at  time  t  (lb/in.2) . 

The  pressure  at  time  t  is  given  by 


(A-5) 


P  -  P0  +  12.0  (cb  -  cex)  FpTk/V 


(A-6) 


where  (cb  -  cex)  is  the  weight  of  gas  in  the  chamber  (lb) ; 

Fp  is  the  impetus  of  the  propellant  (ft-lb/lb) ; 

Tk  is  the  ratio  of  the  isobaric  flame  temperature  (°K)  of 

the  propellant  to  the  isochoric  flame  temperature  (°K) 
of  the  propellant. 


Combining  Equations  A-l  through  A-6, 


12.0 


P  -  Po  + 


z(l  +  ®f)  cp  -  Cd  At  fQ  P  dtj 


FpTk 


Vo  +  z(l  +  0f)  c. 


(A- 7) 


Since  f  «  1  -  z , 


P  -  Po  + 


12.0  ^[-6z2  +  z(l  +  9)  j 


cp  "  ^d  ^t 


£ p  dt) 


FpTk 


VQ  +  [-0z2  +  z(l  +  9)]c. 


(A-8) 


Combining  Equations  A-3  and  A-4, 

w  dz 


dt 


BPn 

av 


(A-9) 
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Integrating , 


(Ar 10) 


Since  z(o)  =  0  and  z(t)  =  z,  and  since  w  and  B  are  constants, 


(A- 11) 


Substituting  Equations  A-9  and  A- 11  into  the  integral  expression 
in  the  numerator  of  Equation  A-8,  and  multiplying  the  numerator  and 
denominator  by  p, 


o 


pV0  +  [-0z*  +  z(l  +  0)]  cp 


pFpTk 

(A-12) 


Equation  A-12  is  Equation  1  in  the  text. 


APPENDIX  B 


ITERATIVE  PROCEDURE  FOR  SOLVING  EQUATION  1 


The  procedure  for  calculating  P  in  Equation  1  is  as  follows: 

1.  At  t  »  0,  P  ■  PQ  and  z  =  0. 

2.  z  ■  z  +  dz. 

3.  Compute  Cb  using  Equation  A-l. 

4.  Compute  V  using  Equation  4- 2'.‘, 

5.  Estimate  a  value  of  P,  making  sure  the  estimated  value  is 
less  than  the  actual  value. 

6.  Compute  r  using  Equation  A-3 . 

7.  Compute  dt  using  Equation  A-4. 

8.  Compute  cex  using  Equation  A- 5. 

9.  Compute  P  using  Equation  A-6. 

10.  If  the  difference  between  the  value  of  P  estimated  in  step 

5  and  the  value  of  P  calculated  in  step  9  is  large,  increment  the  es¬ 
timated  value  of  P  by  some  AP  >  0  and  go  back  to  step  6.  If  this 
difference  is  small,  the  value  of  P  computed  in  9  is  adequate  as  the 
pressure,  P(z  +  dz)  ,  at  the  time  t  +  dt,  where  t  is  the  time  corref* 
sponding  to  the  pressure  P(z)  and  dt  is  the  value  last  computed  in 
step  7.  Proceed  to  step  11. 

11.  If  z  <  1,  go  back  to  step  2;  otherwise,  set  z  «  1  and  pro¬ 
ceed  . 

12.  Slnde  z  ■  1,  the  propellant  is  completely  consumed  and  gas 
is  no  longer  being  generated.  Hence,  the  gas  in  the  chamber  will  be 
expelled  and  the  pressure  will  rapidly  decrease.  The  independent 
variable  is  now  time,  t.  So,  t  ■  t  +  dt,  where  dt  is  designated. 

13.  Set  Cb  ■  cp. 

»  14,  Compute  V  using  Equation  A-2. 

15.  Estimate  a  value  for  P,  making  sure  the  estimated  value  is 
less  than  the  actual  value. 
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16.  Compute  cex  using  Equation  A-5. 

17.  Compute  P  using  Equation  A-6. 

18.  Compare  the  pressures,  as  was  done  in  step  10;  increment  the 
estimated  pressure  and  go  back  to  step  16.  When  good  agreement  is 
obtained,  proceed  to  step  19. 

19.  Increment  t  by  dt,  and  go  back  to  step  15.  Continue  in  this 
manner  until  the  pressures  become  relatively  low;  then  stop. 

As  indicated  above  and  in  the  schematic  flow  chart  (Figure  B-l), 
the  procedure  is  to  assume  the  web  has  regressed  by  a  fraction  dzi 
estimate  what  the  resulting  pressure  would  be,  and  compare  the  cal¬ 
culated  pressure  with  the  estimated  pressure.  This  comparison  is 
made  by  checking  the  truth  of  the  following  inequality: 

| P  -  Pg  |  <  EPS . 

If  the  inequality  is  false,  a  new  estimated  pressure  is  obtained  by 
incrementing  the  present  estimated  pressure  by  EPS.  This  process  is 
continued  until  the  inequality  holds  true. 

The  question  that  naturally  arises  is:  "What  values  of  dz  and 
EPS  are  required  to  produce  a  satisfactory  simulation?"  Apparently, 
smaller  values  of  dz  and  EPS  are  closer  to  the  limit  than  larger 
values.  However,  smaller  values  of  dz  and  EPS  mean  not  only  more 
iterations  in  finding  the  pressure  corresponding  to  a  given  value  of 
dz,  but  also  more  values  of  dz  to  evaluate.  1 

To  determine  the  values  of  dz  and  EPS  that  could  be  used  to 
produce  a  good  simulation  and  yet  not  require  unrealistic  computer 
time,  a  hypothetical  rocket  motor  was  assumed  and  pressure  vs  time 
curves  were  obtained  for  two  cases:  Case  I,  dz  -  0.01,  EPS  «  25  p«*i; 
Case  II,  dz  ■  0.001,  EPS  ■  2.5  psi.  Figure  B-2  illustrates  the  sig¬ 
nificance  of  using  these  sets  of  variables  for  EPS  and  dz.  For  Case 
II,  ten  equally?*  ppacdd  points ,  .whose  ordinates  are  correct  within 
±2.5  pat,  a*$L. determined  on  the  interval  z  ■  0  to  z  “  0.01.  In 
Case  I,  only  one  point  whose  ordinate  is  correct  within  ±  25  psi 
is  determined. 

Figure  3*3  r'shows:  pressure-time  curves  produced  by  the  computer 
in  Case  I  and  Case  II.  It  can  be  seen  that  the  curves  are  remarkably 
close.  Since  Case  II  requires  a  considerable  amount  of  computer  time, 
the  small  over-all  precision  obtained  is  hardly  warranted.  Thus,  the 
values  of  dz  ■  0,01  and  EPS  ■  25  psi  may  be  considered  sufficient  for 
good  simulation. 

The  computation  of  the  pressure  and  thrust  vs  time  curves  is 
terminated  when  the  values  computed  for  the  thrust  coefficient,  Cf, 
become  negative. 
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G0MPTJTE  PAV.PXT, 
CEX,  GAS,  P 


If  -  p«  |  -  eps; 


|p  -  PC|  -  KPS 


PiRMra  B-l.  Schematic  Plow  Chart 
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Figure  B-3.  Pressure  vs  Time  for 

dz  =  0.001,  EPS  *  2.5  psi 

dz  =  0.01,  EPS  =*  25  psi 


APPENDIX  C 


USE  OF  NEWTON-RAPTHSON  METHOD  IN  SOLVING  FOR  Cf 


The  values  of  €,  y,  pa,  and  jie/P  in  the  pair  of  parametric  equa¬ 
tions  (Equations  3  and  4)  are  considered  constant.  Hence,  the  thrust 
coefficient,  Cf,  varies  with  the  pressure,  P.  To  reduce  these  equa¬ 
tions  into  a  form  easier  to  work  with,  the  following  substitutions 
are  made : 


01  -  MU2 

(C-l) 

g2  “  7 

(C-2) 

Go  0  y  -  i 

63  y 

(C-3) 

Pe 

X  “  — 

P 

(C-4) 

Equation  4  then  becomes 

Gi 

(C-5) 

*  XG2  (: 

1  :\Gvm 

which  may  be  rearranged  into 

the  following  fonn 

2G2 

x  - 

x2G2+G3  .  |flj2 

(C-6) 

or 

xm  _  xn  m  1 

(C-7) 

in  which 

•  #  ■  $ 

3±  1 

M  fchT 

(C-8) 

m  ■  2Go  m  2 

y 

n  -  2G2  +  G3  -  |  (C-9) 
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In  the  physical  situation,  n  is  always  greater  than  m,  and  x 
is  always  less  than  1.  We  are  interested  in  determining  the  value 
of  x  that  makes  Equation  C-7  hold  true,  or,  we  are  looking  for  the 
roots  of  f(x)  in  Equation  C-10. 

f(x)  -  xm  -  xn  -  l  (C-10) 

As  a  typical  example,  consider  the  case  when  y  -  1.2  and  €  =  4. 
Using  Equations  C-8,  C-9,  and  C-10,  respectively,  it  is  determined 
that 


t  =  0.0021905869 
m  =  1.6666667 
n  »  1.8333333 

Figure  C-l  is  a  plot  of  y  -  xm  and  y  =  xn,  in  which  m  and  n 
have  the  above  values.  Figure  C-2  shows  a  plot  of  y  =  xm  -  xn,  a  plot 
of  y  ■  t,  and,  also,  a  plot  of  the  difference  of  these  graphs,  viz, 
f(x)  »  x111  -  xn  -  -t.  It  is  evident  from  this  figure  that  f(x)  has 
two  roots.  The  smaller  root  is  desired  since  it  is  physically  mean¬ 
ingful  . 

To  use  the  Newton-Rapthson  (Ref  4)  method,  we  must  have  an  analyt¬ 
ical  expression  for  the  derivative  of  f(x). 

f'(x)  **  m  xm“l  -  n  xn”l  (C-ll) 

Making  sure  that  the  initial  guess  for  x  is  far  to  the  left  of 
the  point  where  f(x)  crosses  the  abscissa,  compute  the  value  for  f(x) 
at  x  ■  x0l  determine  the  value  of  f(x)/f"(x),  add  this  to  x0  and  re¬ 
peat  the  process  using  this  new  value  for  x  as  dictated  by  Equation  C-12. 

Xn+1  "  xn  +  (c-12) 

The  iteration  ceases  when  f(x)  is  zero  or  very  small. 

Table  C-I  shows  the  result  of  this  procedure  using  the  values 
of  y  and  €  cited  above  and  using  an  initial  guess  of  0.001  for  x. 

Only  four  iterations  were  required  in  determining  the  desired  root. 

Having  employed  the  Newton-Rapthson  method  in  solving  for  x, 
and  recalling  that  x  la  really  pP/P,  the  value  for  Cf  may  be  calcu¬ 
lated  directly  by  moans  of  Equation  3. 
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Figure  C-l.  Plot  of  y  =  x1.6667  aniJ  y  „  x1.8333 


Figure  C-2.  Plot  of  y  -  .  ,1.8333  and  y  .  0.002191  and  y  -  x1-66*?  .  ,1.8333  .  0.002191 


TABLE  C-l.  Values  Obtained  in  a  Newton-Rapthson  Iterative  Procedure 
for  solving  x* -6666667  -  xl. 8333333  _  0.0021905869 
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APPENDIX  D 


THE  QUANTITY  cp  6 


The  dimensionless  quantity  6  (ref  2)  in  Equation  2  is  a  combined 
efficiency  factor  and  nozzle  expansion  factor.  The  value  <p  is  an  em¬ 
pirically  derived  constant  and  is  usually  about  0.96  for  good  nozzle  de¬ 
signs.  The  value  6  is  express  as 

6  =  0.5  (1  +  cos  Ct)  (D-l) 

where  Oi  is  the  half  angle  of  the  expansion  cone. 

Figure  D-l  is  a  plot  of  the  combined  factor  <p  6  as  a  function  of 
Q!  for  values  of  <p  equal  to  0.94,  0.96,  and  1.0 
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Figure  D-l.  Combined  Energy  Loss  Correction,  p,  and  Divergence  Angle  Correction, 
vs  Divergence  Angle,  Ot 


APPENDIX  E 


THE  FORTRAN  PROGRAMS 

Rocket  Motor  Performance  Studies  -  Fortran  I 


c  ROCUT  MOTOR  PERFORMANCE 
c  STUDIES 

REAU.  P 

READ*  CP.4tB.F0RN 
RE AO i  FP.OEN.B.BN.G.TK 
REAUt  VOl. •  AT *E* DECT 
READ i  CD. PA 

REAO.  EPS.UUESS.XO.OZ.TI.AIO 
PR  I NT . P 

PRINT. C^.RtS. FORM 
PRINT .FP.OEN.B.BN.to.TK 
PRINT. VOu. AT. E.OELT 
PRINT* CO. PA 

PRlNT.EP3.UUeSS.X0.0Z.TI 'AID 
OR»UZ*«Ed 

ai»«u-i»/a*o 

cxPRiaun/ia-ii 

G24C2.U/CGT1 >)**£»* 

Au»ui*aa>^iE*e> 

AMR2.0/Q 

anru+ui/b 

CUH<M2.0FG*G/<U-1.UI  I RB2 

cu»«mg-ii/u 


1  FXR«X**AM)-(X**AN|-AL 

IFMABSFIFXn-l.0E-TI2.2i3 

s  fpxr«ah*«x**iam-i i i ) 

1  - < AN* ( X*« ( AN-1 ) 1 1 

hrfx/ppx 


XRX-H 
60  TO  I 

2  OEU»CUS*UONF(XI 
IF (UEUI30.31 • SI 
SU  OELRU.U-OEL 

ZlBl.O/EXPPIUEl.1 
60  TO  52 
31  Z1BEXPFI0ELI 
52  Zl«CUR*(l.O-Zll 
Z2RSURTF ( Z1 t 
M«C»X 
CFMR22*Z3 
CFNKCUPA 
TRO.U 
P6R0.0 
PXTRRO.O 
ZAO* 

41  FRA1.U-* 

CBBZ*U.OT(RORM*PR»  »*CP 


vbvoumcb/oeni 

PBRP 

42  P6RPUTEP3 
PAVRIPBAP6I/2.0 

rateau*ipav*rbn> 

DTRtM/RATC 

bxTRPXTRAIPAV*OTI 

CEX«CO*AT*PXT 

6AS4ICB-CCAI 

?g  PR  ( |2.W«'AAh*PP*TK/y  MPA 
I F I ABSF « P-PO ) -EPS  1 43 * *3 • 

43  PRTRAPAT 
CFRCFM- I CFNPP I 
IFICFI4T.4T.4R 

4«  FORCF*P»AT*OCl.T 


42 


TAT+OT 
TGASRTKUUAS 
PRINT  »2»T.P.F0 
IFM-1.0I6S.66.66 
63  ZHZ+OZ 

PURP-GUESS 
GO  TO  61 

66  2R1.0 

PUAP-IaI0*4U£SS) 

rateau.o 

OTATI 

CBFCP 

VRVOL+ICd/UEN) 

PBRP 

Tl  PGRPGTEPS 

PAVRIPB+PG 1/2.0 
PXTRPXTRT ( PA V*uT I 
CEXACD*A  T  *px  T 
gasrcb-cex 

PR ( 12.0RGASUFP6TK/V I ♦PA 
IF! ABSF I P— PG ) — bPS ) 63 1 63 1  7 1 

67  PAVURPXTX/T 
PEXRPAVQAX 
PRINT. PAVG.PEX.X 

CF  AVRCFM-  <  CFi./PAVG  I 
PRINT. CFAV 

FAVGACFAV*PA  VG*AT*l)ELT 

PRINT. FAVO 

edotrco*at«pavu 

PKI.  T.RDUT 

SPIPFAVG/RUOT 

PRINT. SP1 

TiMP»FAva*r 

PRINT. TIMP 

RAVGRB*(PAVa**8N) 

PRINT. KAYO 
END 
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appendix  f 


EXAMPLE  -  COMPARISON  OF  COMPUTER  AND  EXPERIMENTAL  RESULTS 

The  parameters  used  to  simulate  the  operation  of  the  XM15  Rocket 
Catapult  are  presented  in  Table  F-I. 

TABLE  F-I.  Parameters  Describing  the  XM15  Rocket  Motor 


Po 

=  250  lb/in.2 

At 

= 

19.24  in.2 

CP 

-  82  lb 

e 

= 

4.662 

w 

»  0:88  in. . 

tp  6 

a 

.95 

0 

00 

rH 

• 

n 

CD 

a 

.00691 

F  •• 

-  297,700  ft  lb/lb 

EPS 

a 

25  lb/in.2 

P 

-  0.0576  lb/in.3 

P8 

a 

100  lb/in.2 

B 

=0.086  (in./sec) (in.2/lb)n 

xo 

= 

.001 

n 

CM 

• 

i 

dz 

= 

.01 

y 

=  1.23 

Ti 

= 

.005 

Tk 

-  .815 

AID 

= 

8.5 

V  =  2120  in.3 

Figure  F-l  shows  a  comparison  between  the  calculated  pressure  vs 
time  curve  and  a  pressure  vs  time  curve  obtained  from  a  test  fir  ng 
of  the  XM15  rocket  catapult.  Figure  F-2  is  a  plot  of  the  corres¬ 
ponding  calculated  thrust  vs  time  curve. 

Table  F-II  is  a  listing  of  additional  parameters  computed  at  the 
end  of  the  program. 


TABLE  F-II.  Additional  Computed  Parameters 


Pex/P 

m 

.03303 

w 

■  146.4  lb/sec 

Pcx 

- 

38.38  lb/in. 2 

■^sp 

*  214.6  sec 

P 

m 

1101.4  lb/in.2 

*tot 

-  17,590  lb-sec 

Cf 

m 

1.561 

r 

■  1.629  in. /sec 

F 

m 

3,142  lb 

23 


ooo'os 


© 

o 

o 

«k 

o 

CM 

(qt)  xsmmx 


25 
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Figure  F-2.  Computer  Thrust  vs  Time  Curve,  simulating  XM15  Test  Firing 
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1  -  Dir,  NASA  Langley  Research  Center 
Attn:  Mrs.  E.R. Gilman,  Librarian 
Langley  Station 
Hampton,  Virginia,  23365 


National  Aeronautics  and 
Space  Administration 
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1  -  Director,  NASA  Geroge  C. 

Marshall  Space  Flight  Center 
Attn:  Saturn  System  Office 
Huntsville,  Alabama 

1  -  Attn:  Light  and  Medium 
Vehicles  Office 

1  -  Attn:  Propulsion  and  Energy 
Division 

1  -  Attn:  Future  Projects  Office 
1  -  Attn:  Research  Projects  Office 

1  -  Director,  NASA  Manned  Space¬ 

craft  Center 
Attn:  Gemini  Spacecraft 
Projects  Office 
Houston  1,  Texas 

2  -  Attn:  Mr.  Chas.  M.  Grant,  Jr. 

Chief,  Technical  Infortfia- 
tion  Division,  Code  ACI 

2  -  Attn:  Apollo  Spacecraft 
Projects  Office 

1  -  Attn:  Office  of  Asst  Director, 
Research  *  Development 

1  -  Attn:  Spacecraft  Research  Div 

1  -  Attn:  Crew  System  Division 

1  -  Director 

NASA  Goddard  Space  Flight  Center 
Greenbelt,  Maryland 

1  -  Director 

NASA  Jet  Propulsion  Laboratory 
Pasadena,  Calif. 

1  -  Director 

NASA  Lewis  Research  Center 
Cleveland,  Ohio 


U  .S  .Naw 

1  -  Chief,  Bureau  of  Naval  Weapons 
Department  of  the  Navy 
Attn:  RMMO 
Washington  25,  D.  C. 

1  -  Attn:  RMMP 

1  -  Commanding  Officer 

U.S .Naval  Weapons  Laboratory 
Attn:  Mr.  R.  I.  Rossbacher, 

Code  WC 

Dahlgren,  Virginia 

1  -  Commanding  Officer 

U.S .Naval  Ordnance  Test  Station 
Attn:  D.  H.  Strietzel 
Code  4574 
China  Lake,  Calif. 

1  -  Commanding  Officer 

U.S. Naval  Ordnance  Plant 
Attn:  Code  PD- 270 
Macon,  Georgia 

1  -  Commanding  Officer 
Naval  Air  Depot-Crane 
Attn:  Code  4100 
Crane ,  Indiana 

1  -  Director 

Air  Crew  Equipment  Laboratory 
Naval  Air  Engineering  Center 
Attn:  C.  T.  Koochembere, 

Group  Manqger 
Philadelphia  12,  Pa. 

1  -  Commander 

Naval  Air  Development  Center 
Attn:  Technical  Library 
Johnsvtlle,  Pennsylvania 

U.S .Air  Force 

3  -  Commander 

Aeronautical  Systems  Division 
Attn:  ASNRMF-24,  Mr.  R.  Dobbek 
Wright-Pntterson  AFB ,  Ohio 
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3  -  Commander 

Aeronautical  Systems  Division 
Attn:  ASNMES,  C.  K.  Hodell 
Wright-Patterson  AFB,  Ohio 

1  -  Attn:  ASNMEP ,  Mr.  A.  E.  Varble 

1  -  Attn:  ASNMEP,  C.  E.  Grindstaff 

1  -  Attn:  ASZTE,  Capt  J.  J.  Sollars 

1  -  Attn:  ASZWO,  Capt  F.  H.  Briggs 

1  -  Attn:  ASZFO-A 

1  -  Attn:  ASZHPC,  Capt  W.  C.  Knothe 

1  -  Attn:  ASNRH,  Mr.  Geo.  C.  Dickson 

1  -  Attn:  ASBMB,  Maj  N.  P.  Clarke 

1  -  Attn:  ASWMRD,  Mr.  W.  M.  Cheslak 

1  -  Attn:  ASAPRL 

1  -  Attn:ASZB,  Col  C.  A.  Gayle 

1  -  Directorate  of  Aerospace  Safety 
Attn:  Chief,  Life  Sciences  Group 
AFIAS-L 

Norton  Air  Force  Base,  Calif. 

1  -  Commander 

Ogden  Air  Materiel  Area 

Attn:  OONEAA 

H^ll  Air  Force  Base  Utah 

1  -  Commander 

2705th  Airmunition8  Wing  (AMC.) 
Attn:  OOYTC 

Hill  Air  Force  Base,  Utah 

1  -  Commander 

Edwards  Air  Force  Base 
Attn:  DGSM,  6593d  Test  Croup 
Edwards,  Calif. 

1  -  USAF  Special  Weapons  Center 
Attn:  SWOI 

Klrtland  Air  Force  Base,  N .  M. 


Other 

10  -  Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 

1  -  Chemical  Propulsion  Information 
Agency 

Johns  Hopkins  University 
Silver  Spring,  Maryland 

1  -  Scientific  and  Technical 
Information  Facility 
Attn:  NASA  Representative 
(S-AK/DC) 

P.0.  Box  5700 
Bethesda,  Maryland 

1  -  ICRPG  Working  Group  on  Design 
Automation 

The  Johns  Hopkins ;University 
8621  Georgia  Avenue 
Silver  Springs,  Md. 

1  -  Mr.  Boyd  C.  Brinton 

Propulsion  Systems  Analysis  Division 
Thiokol  Chemical  Corporation 
Wasatch  Division 
Bringham  City,  Utah 
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